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Abstract 

We present a method for determining the rms mass fluctuations on 8 h _1 Mpc scale, 
<7g. The method utilizes the rate of evolution of the abundance of rich clusters of galaxies. 
Using the Press-Schechter approximation, we show that the cluster abundance evolution is 
a strong function of <j%: dlogn/dz oc — 1/of; low o"g models evolve exponentially faster than 
high erg models, for a given mass cluster. For example, the number density of Coma-like 
clusters decreases by a factor of ~ 10 3 from z = to z ~ 0.5 for crg=0.5 models, while the 
decrease is only a factor of ~ 5 for cr 8 ~ 1. The strong exponential dependence on o"g arises 
because clusters represent rarer density peaks in low a 8 models. We show that the evolution 
rate at z <^ 1 is insensitive to the density parameter Q or to the exact shape of the power 
spectrum. Cluster evolution therefore provides a powerful constraint on erg. Using available 
cluster data to z ~ 0.8, we find <7g = 0.83 ± 0.15. This amplitude implies a bias parameter 
b ~ o-g 1 = 1.2 ± 0.2, i.e., a nearly unbiased universe with mass approximately tracing light 
on large scales. 

subject headings : galaxies : clusters - galaxies : evolution - galaxies : formation - cos- 
mology : theory - large scale structure of universe 



1. Introduction 



Density fluctuations in the early universe provide the initial seeds for the structures we 
see today. Various observations of large-scale structure constrain the shape of the power 
spectrum of fluctuations on different scales (Fisher et al. 1993, Park et al. 1994, Peacock 
& Dodds 1994, Strauss & Willick 1996, Landy et al. 1996), but the amplitude of the mass 
fluctuation spectrum has been difficult to obtain. (The microwave background anisotropy 
provides so far the only direct amplitude measure on large scales; Smoot et al. 1992, Bunn & 
White 1997). Knowledge of both the shape and amplitude of the spectrum is needed before 
critical cosmological implications can be derived. 
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The amplitude of the mass fluctuations provides an important constraint on the relative 
distribution of mass and light in the universe, which remains one of the central unresolved 
problems in cosmology: does mass follow light on large scales (an unbiased universe) or is 
mass distributed differently than light (biased universe)? The bias parameter (Kaiser 1984) 
is quantified as b = 5 ga i/S m ~ a 8 (gal)/a 8 — l/cr 8 , where 5 = Ap/p is the overdensity (in 
galaxies or mass, respectively), <t 8 is the rms mass fluctuation within a top-hat radius of 
8h _1 Mpc, and a 8 (gal) ~ 1 is the observed optical galaxy fluctuations on that scale (Davis & 
Peebles 1983). A bias of b ~ a 8 ~ 1 implies an unbiased universe, where mass follows light 
on these scales (assuming a scale-independent bias on large scales), while b ~ 2 (a 8 ~ 0.5) 
implies a highly biased universe, with mass distributed more widely than light. 

Galaxy clusters are useful tools in cosmology. Since clusters represent high density 
peaks, their present-day abundance directly reflects the amplitude of density fluctuations 
on the relevant cluster mass scale; the latter corresponds to the mass enclosed within 
R ~ 8h _1 Mpc spheres in the early universe. Observations of the present-day cluster abun- 
dance indeed place one of the strongest constraints on cosmology and the amplitude of mass 
fluctuations : a 8 n - 5 ~ 0.5 (Bahcall & Cen 1992, White et al. 1993, Viana & Liddle 1996, 
Eke et al. 1996, Pen 1997). This relation, while powerful, is degenerate in Q and a 8 ; neither 
parameter can be directly determined from this constraint. Models with Q = 1 and cr 8 ~ 0.5 
are indistinguishable from models with Q ~ 0.25 and a 8 ~ 1. 

In the present paper we investigate a method that can measure, for the first time, the 
amplitude erg, nearly independent of Q. The method is based on the evolution rate of cluster 
abundance to z ~ 1; this evolution breaks the degeneracy that generally exists between 
Q and a 8 and allows the first determination of each of the parameters independently (see 
also Press & Schechter 1974, Peebles et al. 1989, Henry et al. 1992, Eke et al. 1996, 
Oukbir & Blanchard 1997, Carlberg et al. 1997, Bahcall, Fan & Cen, 1997). In our previous 
paper (Bahcall et al. 1997) we presented results from cosmological simulations. In the 
current paper, we develop the theoretical basis that shows that the evolution rate of cluster 
abundance is strongly dependent on cr 8 and insensitive to Q or to the exact shape of the 
power spectrum (to z <; 1). Therefore a 8 can be determined by observations of cluster 
evolution. We present the method in §2, and determine <r 8 in §3. 

2. Dependence of Cluster Evolution on Cosmological Parameters 

The Press-Schechter (1974) formalism approximates the evolution of the comoving num- 
ber density of halos with mass M (mass function) in an initial Gaussian density field assuming 
hierarchical structure formation from linear perturbations: 
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where <r(z, M) is the linear theory rms mass density fluctuation in spheres of mass M at 
redshift z, 5 C ~ 1.68 is the critical density contrast needed for collapse (only weakly dependent 
on fl, Eke et al. 1996, Kitayama & Suto 1996), and p is the mean cosmic density. The mass 
refers to the virial mass of the system. Here 

a (M)D(z,n,A) 
^' M)= D(z = 0,Q,A) ' (2) 

where .D is the linear growth factor (Peebles 1980) that represents the growth of a pertur- 
bation in the linear regime. For Q = 1, D oc (1 + z)^ 1 ; D grows slower for fl < 1. ctq{M) is 
the present rms mass fluctuation within sphere of mass M, 

° 2 o(M) = ^- 2 J dkk 2 P(k)\W R (k)\ 2 , (3) 

where P(k) is the mass density power spectrum and Wn(k) is the window function. If P(k) 
is approximated as a power law P(k) oc k n around 8 h _1 Mpc scale then 

/ M \ ~ a 

a (M) = a 8 (—) oc a 8 M~ a Q a , (4) 

where a — (n + 3)/6, Mg is the mean mass within a sphere of radius 8h _1 Mpc, and M§ oc fl. 
For a CDM-type spectrum, n m — 1 to —2 on scales of ~ 8h _1 Mpc; for mixed (cold + hot) 
dark matter models, MDM, n 2.5 (Klypin et al. 1993). 

We calculate the cluster evolution by numerically integrating eq.(l) (see below). First, 
however, we show how eq.(l) can be approximated, in order to better understand the major 
factors that dominate the evolution. For clusters of Richness ^ 1 (M > 3 x lO 14 h _1 M ), the 
exponential part of eq.(l) dominates the mass function. Therefore, 

*■ (x exp (^U ,5) 



dM ^\2a(z,M) 
Transferring to an integrated mass function, we find approximately (c.f., Pen, 1997) 

\nn(z,> M) oc —(D(z, fl, A)<j (M))~ 2 

oc -a s 2 tt- 2a D- 2 M 2a (6) 

We define the evolution rate £(z), for clusters with mass > M, as : 

S(z,>M) = - dl ° gw f> >M) (7) 

dz 
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For small redshifts, the cluster abundance decreases approximately as n(z) oc 10 £ ( z ^ 1+z \ 
with n(z)/n(0) = I0~ £ ^ z (see below). For masses of rich clusters, eq.(6) yields: 

S(z, > M) oc a 8 2 n- 2a dD ' 2 ^ ,Q ' A) M 2a (8) 

The evolution of the cluster abundance is therefore controlled by four factors : 

1. <r 8 ; the evolution rate depends exponentially on <7g~ 2 ; clusters in low a 8 models evolve 
exponentially faster than in high <r 8 models (for a given Q). This arises since a fixed mass 
cluster corresponds to a higher relative density peak (i.e., rarer fluctuation) in low <r 8 models. 
This is the dominant factor in the cluster evolution at recent times. 

2. dD~ 2 /dz; a linear perturbation grows slower in a low-density universe than in a high- 
density one. For Q = 1, the fluctuations grow as (1 + z) -1 , thus continuing to evolve strongly 
at low z, while a low-density universe freezes out fluctuations at an earlier epoch. High Q 
models thus evolve faster than lower f2 models in recent times. For the same f2, dD~ 2 /dz is 
larger in a A-dominated universe than in an open universe (Peebles 1980), therefore stronger 
evolution is expected in a A-dominated model (see also §3). 

3. VL~ 2a \ this factor arises because a fixed cluster mass M, corresponds to a different 
part of the power spectrum P(k), since k~ l ~ R ~ (M/fi) 1 / 3 is larger for low f2 models. 
According to eq.(4), <t (M) oc £l a . For a given <r 8 , a fixed mass cluster corresponds to a 
smaller a (M) for a lower Q, implying a rarer peak in the density field. Contrary to factor 
#2 above, the contribution from this term results in faster evolution for low Q. As we 
show below , the contributions from these two terms nearly cancel each other for z £ 1. 
This cancellation can be also estimated directly from the approximate relations at low z, 
dD- 2 /dz oc f(n) ~ fi 4 / 7 ~ fi a57 (Lightman & Schechter 1990), and !T 2a oc IT - 54 (for 
a = 0.27, ot n = —1.4), which approximately cancel each other. 

4. M 2a ; this term is a constant for a given cluster mass. The evolution is stronger for 
more massive clusters. 

We illustrate the above results in Figure 1, where the dependence of the evolution rate 
£{z) on Q and a$ is calculated by numerically integrating eq.(l). The results are presented 
for a A = CDM model at redshifts and 1 for clusters with virial mass M > 3 x lO 14 h _1 M 
(corresponding to Abell richness class > 1, Bahcall & Cen 1993). A spectrum with n=-1.4 is 
used, i.e. a = 0.27 (corresponding to a power spectrum shape parameter T = Qh ~ 0.25, Eke 
et at. 1996, Viana & Liddle 1996). At z ~ 0, the cluster evolution rate is mainly determined 
by the amplitude a s and is nearly independent of f2 (or T, see below). For <r 8 = 1, £{z) ~ 1, 
implying that the abundance of richness > 1 clusters decreases by a factor of ~ 3 from z = 
to z = 0.5 (~ l0°- 5£ ( 2 )); for <7 8 = 0.5, £{z) ~ 3.5, and the number density decreases by a 
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factor of ~ 10 2 for the same redshift interval. £{z) changes by only ~ 10 % when fl changes 
from 0.2 to 1; the contributions from the dD~ 2 /dz and fl~ 2a terms described above (eq.8) 
approximately cancel each other. The cluster evolution rate at z < 1, £(z) oc a^ 2 for fixed 
mass clusters, can therefore be used to determine <t 8 . At z ~ 1, the evolution rate depends 
on both fl and <7 8 , but cr 8 is still the dominant factor. 

Figure 2a shows the dependence of the evolution rate on the slope of the spectrum n. 
For a CDM-type power spectrum, n ~ — 1 to -2, the evolution is nearly independent of the 
shape of the spectrum (for any fl; see also White et al. 1993, Eke et al. 1996, Pen 1997). 
The evolution rate £(z) changes by <^ 20% when n increases from -1 to -2. An MDM model 
(n ~ —2.5) evolves only slightly faster than Standard CDM models at this mass. Figure 
2b shows the dependence of the evolution rate on the cluster mass threshold. From eq.(8) 
we see that £(z) oc <j% 2 M 2a . The strong dependence of the cluster evolution on <r 8 and 
the milder dependence on mass threshold (for a ~ 0.3) indicates that determining <r 8 from 
cluster evolution is not sensitive to reasonable uncertainties in the cluster mass; changing 
the cluster mass by a factor of two causes a ^ 20% change in <r 8 (see Fig. 3 below). 

Figure 3 presents the cluster abundance ratio n(z = 0)/n(z = 0.5) as a function of cr 8 
for clusters with virial mass M > 3 x lO 14 h _1 M (Richness ^ 1), and M > 6 x lO 14 h _1 M 
(Richness ^ 2), using an n=-1.4 power spectrum. The dots represent all values of fl from 
0.2 to 1 for each <t 8 . The figure illustrates the strong dependence of the evolution rate on 
<jg, and the weak dependence on fl. The relation can be approximated for all fl values, as: 

n(0) j 10 a55 / CT l (A + 0), lO - 55 ^ 015 /-! (A = 0), M > 3 x 10 14 
n(0.5) ~ j 10 a8 / CT i (A^0), 10 a8n ° 15 ^l (A = 0), M > 6 x 10 14 

These relations are represented by the solid lines in Fig.3. More generally, the evolution of 
any rich cluster can be approximated as 

log 10 n(0)/n(0.5) = 0.8M 2a /a 2 (A ^ 0), 0.8fl 0A5 M 2a /a\ (A = 0) (10) 

where M 6 is the virial mass threshold in units of 6 x lO 14 h _1 M , and a ~ 0.27 for CDM. 
Observations of cluster abundance from z = to z ~ 0.5 can thus determine <t 8 by either 
directly integrating eq.(l) (Fig. 3-4) or by the approximate empirical relations (9-10). 



3. Determination of a s 

In a recent paper (Bahcall et al. 1997), we studied the evolution of the cluster mass 
function using large scale N-body simulations. Several cosmological models were studied : 
two Standard Cold Dark Matter models (SCDM, fl — 1, cr 8 = 1.05 and 0.53); a low-density 
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open CDM (OCDM, fi=0.35, a s = 0.79); a low-density A-dominated CDM(LCDM, fi=0.40, 
cr 8 = 0.80), and two Mixed Dark Matter models (MDM, Q = 1, <r 8 = 0.60 and cx 8 = 0.67). In 
Figure 4 we plot the relation between the cluster abundance ratio n(z = 0)/n(z = 0.5) and 
ag 2 as obtained from the N-body simulations, for clusters with mass Mi. 5 > 6.3 x 10 14 h _1 M Q 
(within a physical radius of 1.5h -1 Mpc, as used in observations), and compare it with the 
expected evolution (§2) of \ogn(z = 0)/n(z = 0.5) oc a 8 2 . The solid line is calculated 
from the Press-Schechter relation (eq.l.) for A = CDM model(for the mean of all JTs, 
from 0.2 to 1), for the same mass clusters (within physical radius of 1.5h _1 Mpc, for proper 
comparison with the simulations and observations; the comoving radius used in the previous 
paper is replaced here with the physical radius to better reflect the final observations). For 
the Press-Schechter approximation, the mass is scaled from virial to 1.5h _1 Mpc radius as 
follows: the overdensity within the virial radius is numerically calculated for each model 
(e.g., Eke et al. 1996, Oukbir & Blanchard 1997), thus yielding the virial radius for a given 
virial mass. The virial mass in eq.(l) (i.e., mass within the virial radius) is then scaled to the 
1.5h _1 Mpc mass using the observed mass profile of M oc r ' 64 (White et al. 1993, Carlberg 
et al. 1997). The results are insensitive to this transformation since the two masses (radii) 
are comparable. The dashed line in Fig. 4 represents the best-fit to the model simulations; 
it agrees well with the Press-Schechter relation. Both Press-Schechter and the simulation 
show that the evolution rate is mostly dependent on a$ 2 , and is insensitive to f2 and to the 
shape of the power spectrum. Figure 4 also shows that, as expected (§2), the LCDM model 
evolves somewhat faster than OCDM, for have similar Q and cr 8 . 

The strong dependence of the cluster evolution rate on <r 8 provides a powerful method for 
constraining a$ by comparison with the observed evolution of high-mass clusters to z ~ 0.5 — 
1. We use the CNOC/EMSS sample of high-redshift clusters (Carlberg et al. 1997, Luppino 
& Gioia 1995) to determine the observed evolution (see, e.g., Bahcall et al. 1997). We find 
only mild evolution of cluster abundance. For massive clusters with M 15 > 6.3 x lO 14 h _1 M 
(within physical radius of 1.5h _1 Mpc) we find a mean best-fit observed logn(z) oc z relation 
of \ogn(z = 0)/n(z = 0.5) = 0.70 ± 0.3(ler). This observed evolution is compared with the 
£{z) — cr$ prediction in Fig. 4 (shaded band). The comparison yields a powerful constraint on 
<t 8 : cr 8 = 0.83 ± 0.15 (1 a) (see also Bahcall et al. 1997). These results are consistent with 
those obtained by Carlberg et al. (1997). Since the cluster evolution rate is not sensitive to 
the shape of the power spectrum, the above determination applies to both CDM and MDM 
models. Low cr 8 (highly biased) models predict an extremely low abundance of high mass 
clusters at z ~ 0.5, inconsistent with the observed cluster abundance at z ~ 0.5. 

The evolution of cluster abundance provides the first method that can directly deter- 
mine the amplitude of the mass fluctuations a 8 independently of Q or the power spectrum. 
Equivalently, it provides the bias parameter on these scales. The data suggests that we 
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live in a nearly unbiased universe, with b ~ a^ 1 ~ 1.2 ± 0.2, and that mass approximately 
traces light on large scales. The determination of erg breaks the degeneracy that exists be- 
tween og and Q from the present-day cluster abundance(§l references; Bahcall et al. 1997). 
Using the <7g — f2 relation determined by Eke et al. (1996), and the present constraint of 
a 8 = 0.83 ± 0.15, we find = 0.3 ± 0.1. 

We thank D. Eisenstein, B. Paczynski, D. Spergel and M. Strauss for helpful discussions. 
This work is supported by NSF grant AST93-15368, ASC93-18185, and NASA grant NAG5- 
2759. XF acknowledges support from an Advisory Council Scholarship. 
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Figure 1. Dependence of the evolution rate €(z) of M v > 3 x 10 14 h 1 M clusters on £1 and 
cr 8 for A = models (at z ~ and z ~ 1). 
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Figure 2. Dependence of the evolution rate £{z) on the slope of the power spectrum n and 
on the cluster mass threshold M. Solid, dashed and dotted lines correspond to Q = 1, 0.5 
and 0.3 flat models, respectively. 




Figure 3. Cluster abundance ratio n(z = 0)/n(z = 0.5) versus <7 8 for clusters with virial mass 
M v > 3 x lO 14 h _1 M (Richness ^ 1), and M v > 6 x lO^h^M© (Richness ^ 2), from the 
Press-Schechter relation eq.(l) (represented by dots for all values of £1 from 0.2 to 1 (bottom 
to top) for each a s ). The lines are the best-fit empirical relations (eq.9, with VI ~ 0.5). 
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Figure 4. Cluster abundance ratio n(z = 0)/n(z = 0.5) versus cr 8 from N-body simulations 
(dots), and Press-Schechter approximation (solid line), for M 15 > 6.3 x 10 14 h _1 M Q clusters. 
The observed abundance ratio (§3) is shown by the shaded region. (The solid line is the 
direct P-S relation (eq.l) for the given Mi. 5 mass threshold, for A = and mean of all iTs; 
the dashed line represents the best-fit relation from the simulations.) 



